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Conventional thermo-statistics address infinite homogeneous systems within the canonical ensem- 
ble. However, some 170 years ago the original motivation of thermodynamics was the description 
of steam engines, i.e. boiling water. Its essential physics is the separation of the gas phase from the 
liquid. Of course, boiling water is inhomogeneous and as such cannot be treated by conventional 
thermo-statistics. Then it is not astonishing, that a phase transition of first order is signaled canon- 
ically by a Yang-Lee singularity. Thus it is only treated correctly by microcanonical Boltzmann- 
Planck statistics. This is elaborated in the present paper. It turns out that the Boltzmann-Planck 
statistics is much richer and gives fundamental insight into statistical mechanics and especially into 
entropy. This can be done to a far extend rigorously and analytically. The deep and essential dif- 
ference between "extensive" and "intensive" control parameters, i.e. microcanonical and canonical 
statistics, is exemplified by rotating, self-gravitating systems. In this paper the necessary appear- 
ance of a convex entropy S(E) and the negative heat capacity at phase separation in small as well 
macroscopic systems independently of the range of the force is pointed out. The appearance of a 
critical end-point for the liquid-gas transition in the p — E or V — E phase diagram can be easily 
explained as well the non-existence of a critical end-point of the solid-liquid transition. 



I. INTRODUCTION 

Since the beginning of thermodynamics in the first half 
of the 19. century its original motivation was the descrip- 
tion of steam engines and the liquid to gas transition of 
water. Here water becomes inhomogeneous and develops 
a separation of the gas phase from the liquid, i.e. water 
boils. 

A little later statistical mechanics was developed by 
Boltzmann[l] to explain the microscopic mechanical ba- 
sis of thermodynamics. Up to now it is generally be- 
lieved that this is given by the Boltzmann-Gibbs canon- 
ical statistics. As traditional canonical statistics works 
only for homogeneous, infinite systems, phase separations 
remain outside of standard Boltzmann-Gibbs thermo- 
statistics, which, consequently, signal phase-transitions 
of first order by Yang-Lee singularities. 

It is amusing that this fact that is essential for the 
original purpose of thermodynamics to describe steam 
engines was never treated completely in the past 150 
years. The system must be somewhat artificially split 
into (still macroscopic and homogeneous) pieces of each 
individual phase [2] . The most interesting configurations 
of two coexisting phases cannot be described by a single 
canonical ensemble. Important inter-phase fluctuations 
remain outside, etc. This is all hidden due to the re- 
striction to homogeneous systems in the thermodynamic 
limit. 

Also the second law can rigorously be formulated only 
microcanonically: Already Clausius [3] distinguished be- 
tween external and internal entropy generating mecha- 
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nisms. The second law is only related to the latter mech- 
anism [4] , the internal entropy generation. Again, canon- 
ical Boltzmann-Gibbs statistics is insensitive to this im- 
portant difference. 

For this purpose, and also to describe small systems 
like fragmenting nuclei or non-extensive ones like macro- 
scopic systems at phase-separation, or even very large, 
self-gravitating, systems, we need a new and deeper def- 
inition of statistical mechanics and as the heart of it: of 
entropy. For this purpose it is crucial to avoid the ther- 
modynamic limit. 

The main aspects of this new thermodynamics were 
described in [5-7] and especially were introduced to the 
chemists community in [8, 9]. I will repeat here only the 
basic arguments. Additionally, I will stress the fact that 
negative heat capacity and convex entropy can be seen 
at proper phase transitions of 1. order, i.e. at phase 
separation, in small as well in macroscopic systems inde- 
pendently whether they have long or short range inter- 
actions. As there was a hot discussion at this conference 
about this point, it seems necessary to repeat the argu- 
ments here. 



II. WHAT IS ENTROPY? 

Entropy, S, is the characteristic entity of thermody- 
namics. Its use distinguishes thermodynamics from all 
other physics; therefore, its proper understanding is es- 
sential. The understanding of entropy is sometimes ob- 
scured by frequent use of the Boltzmann-Gibbs canon- 
ical ensemble, and the thermodynamic limit. Also its 
relationship to the second law is beset with confusion be- 
tween external transfers of entropy d e S and its internal 
production aliS. 

The main source of the confusion is of course the lack 
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of a clear microscopic and mechanical understanding of 
the fundamental quantities of thermodynamics like heat, 
external vs. internal work, temperature, and last not 
least entropy, at the times of Clausius and possibly even 
today. 

Clausius [3] defined a quantity which he first called 
the "value of metamorphosis" , in German " Verwand- 
lungswert" in [3]. Eleven years later he [10] gave it the 
name "entropy" S: 



S a 



r b dE 
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(1) 



where T is the absolute temperature of the body when 
the momentary change is done, and dE is the increment 
(positive resp. negative) of all different forms of energy 
(heat and potential) put into resp. taken out of the sys- 
tem. 

From the observation that heat does not flow from cold 
to hot (see section III, formula 8, however section IV) he 
went on to enunciate the second law as: 



<t f > o, 



(2) 



which Clausius called the uncompensated metamorphosis. 
As will be worked out in section IV the second law as 
presented by eq.(2) remains valid even in cases where 
heat flows from low to higher temperatures. 

Prigogine [4], c.f. [2], quite clearly stated that the vari- 
ation of S with time is determined by two, crucially dif- 
ferent, mechanisms of its changes: the flow of entropy 
d e S to or from the system under consideration; and its 
internal production diS. While the first type of entropy 
change d e S (that effected by exchange of heat d e Q with 
its surroundings) can be positive, negative or zero, the 
second type of entropy change diS is fundamentally re- 
lated to its spontaneous internal evolution ("Verwand- 
lungcn", "metamorphosis" [3]) of the system, and states 
the universal irreversibility of spontaneous transitions. It 
can be only positive in any spontaneous transformation. 

Clausius gives an illuminating example in [3]: When 
an ideal gas suddenly streams under isolating conditions 
from a small vessel with volume V\ into a larger one 
(V2 > Vi), neither its internal energy U, nor its tem- 
perature changes, nor external work done, but its in- 
ternal (Boltzmann-)entropy Si eq.(3) rises, by AS* = 
iVln(V2/Vi) . Only by compressing the gas (e.g. isen- 
tropically) and creating heat AE = E x [(V 2 /Vi) 2/3 - 1] 
(which must be finally drained) it can be brought back 
into its initial state. Then, however, the entropy change 
in the cycle, as expressed by integral (2), is positive 
(= TV In (V2/V1)). This is also a clear example for a mi- 
crocanonical situation where the entropy change by an 
irreversible metamorphosis of the system is absolutely 
internal. It occurs during the first part of the cycle, 
the expansion, where there is no heat exchange with the 
environment, and consequently no contribution to the 
integral(2). The construction by eq.(2) is correct though 



artificial. After completing the cycle the Boltzmann- 
entropy of the gas is of course the same as initially. All 
this will become much more clear by Boltzmann's micro- 
scopic definition of entropy, which will moreover clarify 
its real statistical nature: 

Boltzmann[l] later defined the entropy of an isolated 
system (for which the energy exchange with the environ- 
ment d e Q = 0) in terms of the sum of possible configura- 
tions, W, which the system can assume consistent with 
its constraints of given energy and volume: 



S=k*lnW 



(3) 



as written on Boltzmann's tomb-stone, with 



W{E, N, V) 
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(4) 

in semi-classical approximation. E is the total energy, N 
is the number of particles and V the volume. Or, more 
appropriate for a finite quantum-mechanical system: 

rrr/r, , f xr\ \ ^ all cigenstates n of H with given N,V, /r\ 

W{±b,l\, V ) — ^ and E < E n < E + e W 

and eo ~ the macroscopic energy resolution. This is still 
up to day the deepest, most fundamental, and most sim- 
ple definition of entropy. There is no need of the ther- 
modynamic limit, no need of concavity, extensivity and 
homogeneity. In its semi-classical approximation, eq.(4), 
W(E, N, V, ■ ■ ■ ) simply measures the area of the sub- 
manifold of points in the 67V-dimensional phase-space (T- 
space) with prescribed energy E, particle number N, vol- 
ume V, and some other time invariant constraints which 
are here suppressed for simplicity. Because it was Planck 
who coined it in this mathematical form, I will call it the 
Boltzmann-Planck principle. It is further important to 
notice that S(E, N, V) is everywhere analytical in E [11]. 
In the microcanonical ensemble are no "jumps" or multi- 
valuedness in S(E), independently of whether there are 
phase transitions or not, in clear contrast to the canoni- 
cal S(T, N,V). A fact which underlines the fundamental 
role of microcanonical statistics. 

The Boltzmann-Planck formula has a simple but deep 
physical interpretation: W or S measure our ignorance 
about the complete set of initial values for all 6N mi- 
croscopic degrees of freedom which are needed to specify 
the V-body system unambiguously [12]. To have com- 
plete knowledge of the system we would need to know 
(within its scmiclassical approximation (4)) the initial 
positions and velocities of all N particles in the system, 
which means we would need to know a total of 67V values. 
Then W would be equal to one and the entropy, S, would 
be zero. However, we usually only know the value of a 
few parameters that change slowly with time, such as the 
energy, number of particles, volume and so on. We gener- 
ally know very little about the positions and velocities of 
the particles. The manifold of all these points in the 67V- 
dim. phase space, consistent with the given macroscopic 
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constraints of E, N, V, ■ ■ ■ , is the microcanonical ensem- 
ble, which has a well-defined geometrical size W and, by 
equation (3), a non- vanishing entropy, S(E, N, V, ■ ■ ■ ). 
The dependence of S(E, N, V, ■ ■ ■ ) on its arguments de- 
termines completely thermostatics and equilibrium ther- 
modynamics. 

Clearly, Hamiltonian (Liouvillean) dynamics of the 
system cannot create the missing information about the 
initial values - i.e. the entropy S(E, N,V,---) cannot de- 
crease. As has been further worked out in [13, 14] the 
inherent finite resolution of the macroscopic description 
implies an increase of W or S with time when an external 
constraint is relaxed. Such is a statement of the second 
law of thermodynamics, which requires that the internal 
production of entropy be positive for every spontaneous 
process. Analysis of the consequences of the second law 
by the microcanonical ensemble is appropriate because, 
in an isolated system (which is the one relevant for the 
microcanonical ensemble), the changes in total entropy 
must represent the internal production of entropy, see 
above, and there are no additional uncontrolled fluctuat- 
ing energy exchanges with the environment. 



III. THE ZERO'TH LAW IN CONVENTIONAL 
EXTENSIVE THERMODYNAMICS 

In conventional (extensive) thermodynamics thermal 
equilibrium of two systems (1 & 2) is established by 
bringing them into thermal contact which allows free en- 
ergy exchange. Equilibrium is established when the total 
entropy 



S 1+2 (E,E 1 ) = S 1 (E 1 ) + S 2 (E-E 1 ) 



(6) 



is maximal. Under an energy flux AE 2 ^i from 2 — > 1 
the total entropy changes to lowest order in AE by 



AS 1+2 \ E = (T 2 -T 1 )AE 2 ^ 1 . 



(7) 



Consequently, a maximum of S to tai(E, E\)\e > S\ +2 will 
be approached when 



sign(AS, 



total j 



sign(T 2 - T 1 )sign(AS 2 ^ 1 ) > 0. (8) 



From here Clausius' first formulation of the Second Law 
follows: "Heat always flows from hot to cold". Essential 
for this conclusion is the additivity of S under the split 
(eq.6). There are no correlations, which are destroyed 
when an extensive system is split. Temperature is an 
appropriate control parameter for extensive systems. 

It is further easy to see that the heat capacity of an 
extensive system with S(E,N) = 2S(E/2,N/2) is neces- 
sarily positive 



C V {E) =dE/0T=- 



(dS/dE) 2 
d 2 S/dE 2 



> 



(9) 



The combination two pieces of N/2 particles each, but 
with different energy per particle, one at e a = e 2 — Ae/2 



and a second at e b = e 2 + Ae/2, must lead to S(E 2 , N) > 
S(E a /2,N/2) + S(E b /2,N/2), the simple algebraic sum 
of the individual entropies because by combining the two 
pieces one normally looses information. This, however, 
is equal to [S(E a ,N) + S{E b ,N)]/2, thus S(E 2 ,N) > 
[S(E a ,N) + S(E b ,N)]/2. I.e. the entropy S(E,N) of 
an extensive system is necessarily concave, d 2 S/dE 2 < 
and eq. 9 follows. In the next section we will see that 
therefore extensive systems cannot have phase separation, 
the characteristic signal of transition of first order. 



IV. NO PHASE SEPARATION WITHOUT A 
CONVEX, NON-EXTENSIVE S(E) 

At phase separation the weight e s ^~ E / T of the con- 
figurations with energy E in the definition of the canon- 
ical partition sum 



Z(T) 
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e S(E)-E/T dE 



(10) 



becomes bimodal, at the transition temperature it has 
two peaks, the liquid and the gas configurations which 
are separated in energy by the latent heat. Consequently 
S(E) must be convex (d 2 S/dE 2 > 0, like y = x 2 ) and 
the weight in (10) has a minimum at E min between the 
two pure phases. Of course, the minimum can only be 
seen in the microcanonical ensemble where the energy 
is controlled and its fluctuations forbidden. Otherwise, 
the system would fluctuate between the two pure phases 
by an, for macroscopic systems even macroscopic, energy 
AE <~ Ei a t of the order of the latent heat. The heat 
capacity is 



C v (E min )=dE/dT = 



(dS/dEf 
d 2 S/dE 2 



< 0. 



(11) 



I.e. the convexity of S{E) and the negative heat ca- 
pacity are the generic and necessary signals of phase- 
separation[5]. It is amusing that this fact, which is essen- 
tial for the original purpose of Thermodynamics to de- 
scribe steam engines, seems never been really recognized 
in the past 150 years. However, such macroscopic energy 
fluctuations and the resulting negative specific heat are 
already early discussed in high-energy physics by Carlitz 
[15]. 

The existence of the negative heat capacity at phase 
separation has a surprising but fundamental conse- 
quence: Combining two equal systems with negative heat 
capacity and different energy per particle, they will re- 
lax with a flow of energy from the lower to the higher 
temperature! This is consistent with the naive picture 
of an energy equilibration. Thus Clausius' "energy flows 
always from hot to cold", i.e. the dominant control-role 
of the temperature in thermo- statistics as emphasized by 
Hertz [16] is violated. Of course this shows quite clearly 
that unlike to extensive thermodynamics the temperature 
is not the appropriate control parameter in non-extensive 
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situations like e.g. at phase separations, nuclear frag- 
mentation, or stellar systems. [17] 

By the same reason the well known paradox of Antonov 
in astro-physics due to the occurrence of negative heat 
capacities must be reconsidered: By using standard ar- 
guments from extensive thermodynamics Lyndcn-Bell 
[18] claims that a system a with negative heat capac- 
ity C a < in gravitational contact with another b with 
positive heat capacity C& > will be unstable: If initially 
T a > T b the hotter system a transfers energy to the colder 
b and by this both become even hotter! If > — C a , 
T a rises faster than T\> and if the heat capacities don't 
change, this will go for ever. This is Lyndcn-Bclls gravo- 
thermal catastrophe. This is wrong because just the op- 
posite happens, the hotter a may even absorb energy from 
the colder b and both systems come to equilibrium at the 
same intermediate temperature c.f. [17, 19]. Negative 
heat can only occur in the microcanonical ensemble. 

As phase separation exists also in the thermodynamic 
limit, by the same arguments as above the curvature of 
S(E) remains convex, d 2 S/{dE) 2 > 0. Conseguently, 
the negative heat capacity should also be seen in ordinary 
macroscopic systems studied in chemistry! see section 
VIII. 

Searching for example in Guggenhcims book [2] one 
finds some cryptic notes in §3 that the heat capacity of 
steam at saturation is negative. No notice that this is 
the generic effect at any phase separation! Therefore let 
me recapitulate in the next section how chemists treat 
phase separation of macroscopic systems and then point 
out why this does not work in non-extensive systems like 
fragmenting nuclei, at phase separation in normal macro- 
scopic systems, or large astronomical systems. 

V. MACROSCOPIC SYSTEMS IN CHEMISTRY 

Systems studied in chemical thermodynamics consist 
of several homogeneous macroscopic phases ai,ct2,--- 
cf.[2]. Their mutual equilibrium must be explicitly con- 
structed from outside. 

Each of these phases are assumed to be homogeneous 
and macroscopic (in the "thermodynamic limit" (N a — > 
oo | Pa= const))- There is no common canonical ensemble 
for the entire system of the coexisting phases. Only 
the canonical ensemble of each phase separately becomes 
equivalent in the limit to its microcanonical counterpart. 

The canonical partition sum of each phase a is defined 
as the Laplace transform of the underlying microcanoni- 
cal sum of states W(E) a = e s «^ [20, 21] 

Z a (T a ) = / e s ^~ E ^dE. (12) 
Jo 

The mean canonical energy is 

< E a (T a ) >= T 2 a d\nZ a {T a )/dT a . (13) 

In chemical situations proper the assumption of homoge- 
neous macroscopic individual phases is of course accept- 



able. In the thermodynamic limit (N a — > oo| Pcv=ccmst ) of 
a homogeneous phase a, the canonical energy 
<E a (T a )> becomes identical to the microcanonical en- 
ergy E a when the temperature is determined by T^ 1 = 
8S a (E,V a )/dE a . The relative width of the canonical 
energy is 

A^ g ^%H M« ' (14) 

The heat capacity at constant volume is 

< E'i. >t„ — < E a ~>% 
C a \v a = ->0- (15) 

Only in the thermodynamic limit (N a — > oo| Pa=cons t) 
does the relative energy uncertainty AE a — > 0, and the 
canonical and the microcanonical ensembles for each ho- 
mogeneous phase (a) become equivalent. This equiva- 
lence is the only justification of the canonical ensemble 
controlled by intensive temperature T, or chemical poten- 
tial n, or pressure P. I do not know of any microscopic 
foundation of the canonical ensemble and intensive con- 
trol parameters apart from the limit. 

The positiveness of any canonical Cy (T) or Cp(T) (15) 
is of course the reason why the inhomogeneous system of 
several coexisting phases (ai&o^) with an overall nega- 
tive heat capacity cannot be described by a single com- 
mon canonical distribution [5, 22]. 

VI. NEW KIND OF PHASES WELL SEEN IN 
HOT NUCLEI OR MULTI-STAR SYSTEMS. 

The new lesson to be learned is that if one defines the 
phases by individual peaks [35] in e s{E) ~ E/T in (10), then 
there exist also inhomogeneous phases like in fragmented 
nuclei or stellar systems. The general concept of thermo- 
statistics becomes enormously widened. 

Now, certainly neither the phase of the whole multi- 
fragmented nucleus nor the individual fragments them- 
selves can be considered as macroscopic homogeneous 
phases in the sense of chemical thermodynamics (ChTh). 
Consequently, (ChTh) cannot and should not be applied 
to fragmenting nuclei and the microcanonical description 
is ultimately demanded. This becomes explicitly clear 
by the fact that the configurations of a multi-fragmented 
nucleus have a negative heat capacity at constant vol- 
ume Cy and also at constant pressure Cp (if at all a 
pressure can be associated to nuclear fragmentation [5]). 
Meanwhile experimental evidences of negative heat ca- 
pacities have accumulated: Nuclear fragmentation e.g. 
[23], atomic clusters e.g. [24], astrophysics e.g. [25], 
conventional macroscopic systems at phase separation 
e.g. [2]. 

The existence of well defined peaks (i.e. phases as de- 
fined above) in the event distribution of nuclear frag- 
mentation data is demonstrated very nicely in [26] from 
various points of view. A lot more physics about the 
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mechanism of phase transitions can be learned from such 
studies. 



VII. APPLICATION IN ASTROPHYSICS 

The necessity of using "extensive" instead of "inten- 
sive" control parameter is explicit in astrophysical prob- 
lems. E.g.: for the description of rotating stars one con- 
ventionally works at a given temperature and fixed an- 
gular velocity ft c.f. [27]. Of course in reality there is 
neither a heat bath nor a rotating disk. Moreover, the 
latter scenario is fundamentally wrong as at the periph- 
ery of the disk the rotational velocity may even become 
larger than velocity of light. Non-extensive systems like 
astro-physical ones do not allow a "field-theoretical" de- 
scription controlled by intensive fields ! 

E.g. configurations with a maximum of random energy 

en 2 

E r andom E ^ ^pot (16) 

and consequently with the largest entropy are the ones 
with smallest moment of inertia 6 compact single stars. 
Just the opposite happens when the angular-momentum 
L and not the angular velocity O are fixed: 

L 2 

Erandom E ^pot- (1^) 

Then configurations with large moment of inertia are 
maximizing the phase space and the entropy. I.e. even- 
tually double or multi stars are produced, as observed in 
reality. 

In figure 1 one clearly sees the rich and realistic mi- 
crocanonical phase-diagram of a rotating gravitating sys- 
tem controlled by the "extensive" parameters energy and 
angular-momentum. [28] 

VIII. NEGATIVE HEAT CAPACITY AT 
PHASE-SEPARATION CAN ALSO BE SEEN IN 
MACROSCOPIC SYSTEMS INDEPENDENTLY 

OF THE RANGE OF THE INTERACTION. 

The "convex intruder" in S(E) with the depth 
AS sur f(E min ) has a direct physical significance: Its 
depth is the surface entropy due to constraints by the ex- 
istence of the inter-phase boundary between the droplets 
of the condensed phase and the gas phase and the cor- 
responding correlation. AS sur f(E m i n ) is directly related 
to the surface tension per surface atom (with number 
N sur f) of the droplets. 

Irp AS sur f(E m i n ) , 1Q s 

&surf/J-tr = 77 (J-8J 

1* surf 

In my paper together with M.Madjet [29] we have com- 
pared the values of AS sur f(E min ) calculated by Monte- 
Carlo using a realistic short range interaction with the 
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FIG. 1: Phase diagram of rotating self-gravitating systems 
in the energy-angular-momentum (E, L)-plane. DC: region 
of double-stars, G: gas phase, SC: single stars. In the mixed 
region one finds various exotic configurations like ring-systems 
in coexistence with gas, double stars or single stars. In this 
region of phase-separation the heat capacity is negative and 
the entropy is convex. The dashed lines E — L — 1 (left) 
and E = L (right) delimit the region where calculations were 
carried out. 



values of the surface tension of the corresponding macro- 
scopic system. Table (I) shows the scaling behavior of 
AS sur f(E min ) with the size N of the system. 

TABLE I: Parameters of the liquid-gas transition of small 
sodium clusters (in the MMMC-calculation [5, 29] several 
clusters coexist) in comparison with the bulk for a rising 
number N of atoms, N sur f is the average number of sur- 
face atoms (estimated here as Yl ^Huster) OI au clusters with 
Ni > 2 together. a/T tr = As sur f * N/N sur f corresponds to 
the surface tension. Its bulk value is adjusted to agree with 
the experimental values of the a s parameter which we used 
in the liquid-drop formula for the binding energies of small 
clusters, c.f. Brechignac et al. [30], and which are used in this 
calculation [5] for the individual clusters. 





N 


200 


1000 


3000 


bulk 




Ttr [K] 


940 


990 


1095 


1156 




qiat [eV] 


0.82 


0.91 


0.94 


0.923 


Na 


Sboil 


10.1 


10.7 


9.9 


9.267 






0.55 


0.56 


0.44 






■^V surf 


39.94 


98.53 


186.6 


oo 




O/Ttr 


2.75 


5.68 


7.07 


7.41 



Roughly ASsurf(E) oc N 2 / 3 and one may argue that 
this will vanish compared to the ordinary leading volume 
term S vo i(E) oc N. However, this is not so as S vo i(E) at 
energies inside the phase-separation region (the convex 
intruder) is the concave hull of S(E) (its slope gives the 
Maxwell construction of the caloric curve T{E)). It is a 
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FIG. 2: Microcanonical Monte-Carlo (MMMC) [5, 29] simu- 
lation of the entropy s(e) per atom (e in eV per atom) of a 
system of N = 1000 sodium atoms at an external pressure 
of 1 atm. At the energy e < ei the system is in the pure 
liquid phase and at e > ez in the pure gas phase, of course 
with fluctuations. The latent heat per atom is qiat = ez — ei. 
Attention: the curve s(e) is artificially sheared by subtracting 
a linear function 25 + e * 11.5 in order to make the convex 
intruder visible. s(e) is always a steep monotonia rising func- 
tion. We clearly see the global concave (downwards bending) 
nature of s(e) and its convex intruder. Its depth is the entropy 
loss due to additional correlations by the interfaces. It scales 
oc A -1 ' 3 . From this one can calculate the surface tension 
per surface atom a sur f /Ttr = As sur f * N/N sur f. The double 
tangent (Gibbs construction) is the concave hull of s(e). Its 
derivative gives the Maxwell line in the caloric curve e(T) at 
Ttr- In the thermodynamic limit the intruder would disappear 
and s(e) would approach the double tangent from below, not 
of course S(E), which remains deeply convex: The probability 
of configurations with phase-separations is suppressed by the 
(infinitesimal small) factor e~ N 1 relative to the pure phases 
and the distribution remains strictly bimodal in the canonical 
ensemble in which the region ei < e < ez of phase separation 
gets lost. 



IX. THE MICROSCOPIC ORIGIN OF THE 
CONVEXITIES OF S(E) AND OF 
PHASE-SEPARATION 

Many applications of microcanonical thermodynam- 
ics to realistic examples of hot nuclei, atomic clus- 
ters, and rotating astrophysical systems have been pre- 
sented during the past twenty years which demon- 
strate convex intruders in the microcanonical entropy 
and, consequently, negative heat capacities. Such 
are reviewed in the publication list on the web site 
http:/ /www. hmi.de/people/gross/ and elsewhere [3 1-33]. 
Here we shall illuminate the general microscopic mecha- 
nism leading to the appearance of a convex intruder in 
S(E, V, N, , • • • ) as far as possible by rigorous and ana- 
lytical methods. This is the generic signal of phase tran- 
sitions of first order and of phase-separation within the 
microcanonical ensemble. Assume the system is classical 
and obeys the Hamiltonian: 



N 



^ 2m 11 n 
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<S> mt [{r}} ■= 5>(r,- r,-) 

i<3 

In this case the system is controlled by energy and vol- 
ume. For what follows it is important to stress that 
$ m * [{ r }] is independent of the energy as control pa- 
rameter. The topology of its attractive pockets is also 
independent of the available volume V. 



A. Liquid-gas transition 

The microcanonical sum of states or partition sum is: 

1 



W(E,N, V) 



N\(2irK) 3N 

N 



(21) 



( d 3N 7 f d 3N P, e 6(E-Y — - * mt [{ r }]) 
Jv N J 2m * 



straight line and its curvature d S vo i/dE 2 = 0. Conse- 
quently for large N 



d 2 S/dE 2 ~ d 2 S vol /8E 2 + d 2 AS surf /dE 2 + ■■■ 

x d 2 AS surf /dE 2 (19) 



V N e Q {E-E^ N - 2 )/ 2 Yf 1 mT 



I 

Jv N 



N\T{m/2){2TTh 2 f N / 2 
d 3N r ( E- $ mt [{?}] 



_> \ (3JV-2)/2 



yN 



E — Eq 



(22) 



and the depth of the intruder AS sur f(E min ) — a/T tr * 

Nsurf ~ iV 2 / 3 goes to infinity in the thermodynamic Wid- gas (E — E , N, V) x W con f(E — E a , N, V) 

limit. Of course, the ubiquitous phenomena of phase sep- 
aration exist only by this reason. It determines the (neg- 
ative) heat capacity as in eq.(ll). The physical (quite = e l s ^- ga s+s aonf ] ^3) 
surprising) consequences are discussed in [7, 9] . 
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N 3/2 



W id . gas (E,N,V) = 



V N eoE (3N-2)/2Yl* 

N\T(ZN/2){2TTh 2 y N / 2 



(24) 



W conf (E-E Q ,N,V) 



d 3N i 
V N 



Q(E- $ mt [{r}]) 



$ mt [{ r }] - £o 
E — Eq 



(3JV-2)/2 



(25) 



V is the spatial volume; £7 = min <f> m * [{ r*}] is the en- 
ergy of the ground-state of the system. The separation of 
W(E, N, V) into Wid- ga s an d W con f is the microcanon- 
ical analogue of the split of the canonical partition sum 
into a kinetic part and a configuration part: 



V N ( mT\ 



3N/2 



Nl { 2-Kh 2 J 



d 3N r 

yN 



Ml 



(26) 



In the thermodynamic limit, the order parameter of 
the (homogeneous) liquid-gas transition is the density. 
The transition is linked to a condensation of the system 
towards a larger density controlled by pressure. For a fi- 
nite system, we expect analogous behavior. However, for 
a closed finite system, which is allowed to become inho- 
mogeneous at phase separation, this is controlled by the 
available system volume V and not by intensive density 
or pressure. At low energies, the N particles conden- 
sate into a droplet with much smaller volume V ,n "C V. 
3(N — 1) internal coordinates arc limited to Vq,n- Only 
the center of mass of the droplet can move freely in V 
(remember we did not fix the center-of-mass in equation 
eq.(21)). The system does not fill the 37V-configuration 

1/3 

space Vjv- Only a stripe with width V 0N in 3(N — 1) 
dimensions of the total 3iV-dim space is populated. The 
system is non-homogeneous even though it is cquilibrizcd 
and, at low energies, internally in the single liquid phase; 
and it is not characterized by an intensive homogeneous 
density. In fact, W con f(E — Eq,N, V) can be written as: 



W conf (E-E ,N,V) = 

iJV def 



V(E,N) 



V 



N 



< 1 (27) 



[V(E,N)Y 

( d 3N r 0{E- $" l *[{7}]) 
Jv N 



1 $ mt [{r}}-E 

E — Eq 



(3JV-2)/2 



S conf (E-E ,N,V) = Nln 



V{E,N) 



V 



(28) 



< (29) 



The first factor 0(£-$ mt [{7}]) in eq(28) eliminates the 
energetically forbidden regions. Only the potential holes 
(clusters) in the 3iV-dim potential surface $ m *[{r}] < E 
remain. Their volume V N (E,N) < V N is the accessible 



part of the 3iV-dim-spatial volume where $ mt [{r}] < E. 
I.e. V N (E,N) is the total 37V-dim. eigen-volume of the 
condensate (droplets), with N particles at the given en- 
ergy, summed over all possible partitions, clusterings, in 
3iV-configuration space. The relative volume fraction of 
each partition compared with V N (E, N) gives its relative 
probability. V N (E,N) has the limiting values: 



[V(E,N)} N = ) jr N _, 



V N for E in the gas phase 

Von"-^ for E = E 



W conf (E -Eq,N, V) and S conf (E — Eq, N, V) have the 
limiting values: 

W conf (E - Eq) < 1,=> S conf {E - Eq, N) < 

f 1 E > $ m * 



E — > Eq 



(30) 



S C onf(E — Eq) 



ln{[ty 



1JV- 



E > $ m * 

- 1 } < E^ Eq 

(31) 

All physical details are contained in W con f(E — 
Eq,N,V) alias Nln[V(E, N)], c.f. eqs. (27-31): If the en- 
ergy is high the detailed structure of $ mt [{ r }] is unim- 
portant W con f ~ 1, S con f w 0. The system behaves 
like an ideal gas and fills the volume V. At sufficiently 
low energies only the minimum of <I> m * [{ r }] is explored 
by W con f {E — Eq , N, V) . The system is in a condensed 
phase, a single liquid drop, which moves freely inside the 
empty larger volume V, the 3(N — 1) internal degrees of 
freedom are trapped inside the reduced volume Vqn <C V. 

One can guess the general form of N ln[V(E, N)]: Near 
the groundstate E > E it must be flat s=s (N—l) ln[Vbiv] + 
ln[y — Vojv] because the liquid drop has some eigen- 
volume Vqn in which each particle can move (liquid). 
With rising energy ln[V(E,N)] rises up to the point 
(Etrans) where it is possible that the drop fissions into 
two. Here an additional new configuration opens in 37V- 
dim configuration space: Either one particle evaporates 
from the cluster and explores the external volume V , or 
the droplet fissions into two droplets and the two CM co- 
ordinates explore the larger V. This gives a sudden jump 
in Sconf 

(E) by something like ~ ln{ y^'" } and sim- 
ilar jump upwards in the second case. 

Later further such "jumps" may follow. Each of these 
"jumps" induce a convex upwards bending of the total 
entropy S(E) (eq.23). Each is connected to a bifurca- 
tion and bimodality of e s ( E )- E / T an d the phenomenon 
of phase- separation. 

In the conventional canonical picture for a large num- 
ber of particles this is forbidden and hidden behind the 
familiar Yang-Lee singularity of the liquid to gas phase 
transition. In the microcanonical ensemble this is ana- 
logue to the phenomenon of multi-fragmentation in nu- 
clear systems [5, 34]. This, in contrast to the mathemat- 
ical Yang-Lee theorem, physical microscopic explanation 
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of the liquid to gas phase transition sheds sharp light on 
the physical origin of the transition, the sudden change 
in the inhomogeneous population of the 3iV-dim. config- 
uration space. 

B. Solid-liquid transition 

In contrast to the liquid phase, in the crystal phase a 
molecule can only move locally within its lattice cage of 
the size d 3 , instead freely in the whole volume Von of 
the liquid condensate. These are deep localized holes in 
the potential surfaces of <&'"*[{ r }]. I.e. in equation (31) 
instead we have S con f — ► ln{[-^-] The convexity 
of S con f (E) is not controlled by the available volume V 
of the system. 

C. Conclusion 

Only by treating the many-body system microcanoni- 
cally the ~ 170 years old challenge of thermodynamics is 
microscopically solved: 

The essential differences between the gas, the liquid, 
and solid phase are the following: Whereas the gas oc- 
cupies the whole container, the liquid is confined to a 
definite condensate volume, however this may have any 
shape. It is separated from the gas by a surface. The solid 
is also confined to definite volume but in contrast to the 
liquid its surface has also a definite shape. Of course this 
is the standard experimental identification of the phase 
transition and not any singularity (Yang-Lee). However, 



these differences cannot be seen in the canonical ensem- 
ble. As demonstrated in section VIII, microcanonically 
this can well be seen even in macroscopic systems with 
any kind of interaction, short- and long range. 

The gas- liquid transition is linked to the transition 
from uniform filling of the container volume V by the gas 
to the smaller eigen- volume of the system Vo in its con- 
densed phase where the system is inhomogeneous (some 
liquid drops inside the larger empty volume V). First 
3(N — 1), later at higher energies less and less degrees of 
freedom condensate into the drop. First three, then more 
and more degrees of freedom (center-of-mass-coordinates 
of the drops) explore the larger container volume V lead- 
ing to upwards jumps (convexities) of S con f(E). The 
volume of the container controls how close one is to 
the critical end-point of the transition, where phase- 
separation disappears. Towards the critical end-point, 
i.e. with smaller V, the jumps ln[V — Vo] — In [Vo] be- 
come smaller and smaller. At the surface of a drop 
<f> m * > E = min$ mt , i.e. the surface gives a nega- 
tive contribution to S con f in equation (28) and to S at 
energies E>E Q . 

In the case of the solid-liquid transition, however, 
the external volume, V, of the container confines only 
the center-of-mass position of the crystal, resp., the 
droplet. The entropy jumps during melting by AS con f oc 
ln[Vojv] — lnd 3 . These jumps are not controlled by the 
external volume V resp. the external pressure. There is 
no critical end-point of the solid-liquid transition-line in 
the V — E or P — E phase diagram. 

Discussions with S.Ruffo are gratefully acknowledged. 
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